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INTRODUCTION 


ERDA-NASA has planned a graduated wind turbine project, with turbines 
ranging from 100-kw to 1500-kw [1]. The project currently has a 100-kw wind 
turbine at the NASA Plum Brook Station near Sandusky, Ohio, which became 
operational in 1975. This wind turbine, designated the Mod-0, is used as a 
practical case study in this paper (Fig. 1). 

The Mod-0 is a horizontal -axis, propeller type machine. The rotor is down- 
wind of the tower and rotates at a constant speed of 40 rpm. The alternator is 
a 125 kva, 3 phase, 60 hz, 1800 rpm, 480 volt, Y-connected synchronous machine. 
Figure 2 shows details of the wind turbine drive train assembly, and yaw system. 

The control of a system which generates power from as unsteady an input as 
the wind presents a formidable problem. Even with constant wind velocity, the 
input at a wind turbine's blades varies due to the cyclical obstruction of the 
wind by the turbine's supporting tower. Control of the pitch of the turbine's 
blades relative to their plane of rotation is perhaps the most promising method 
to stabilize the effects of such an input upon the mechanical stresses of the 
turbine and the electromagnetic transients of the generator. 

The first purpose of this study is to analyze the current control method for 
wind turbine generators, i.e. blade pitch angle control and a slip clutch [2]. 
Secondly, an investigation is conducted via computer simulation of the effects 
of techniques proposed in this paper upon the dynamic stability of the Mod-0. 

SYSTEM EQUATIONS 

Figure 3 is a schematic diagram of the system under study. Speed control 
is used only for the starting up and synchronization operations. The power 
control system Is used for stabilizing the output under load conditions. The 
blade pitch control system is shown in Fig. 4. 



Synchronous Machine Equations 

Park's equation [3], with appropriate modifications, [4,5], are used. 

It is reasonable to assume that the transformer emfs are 

negligible as compared to the speed voltages Also, the synchronous 

generator frequency may be assumed to be the same as the system frequency, aj Q , 
as long as synchronism is retained. 


^ ■ < V r e ,l> /T do 

(1) 

K. !a2 x d- x 3 
” T do V x d 

(2) 

del 

dt Vqo 

(3) 

(Re^xV) - (RV.+x"Vj 
. _ a d q q v a d q q 

“ *1 * x d x q 

(4) 

, _ <*,'S- X ?5) - (R «V*dV 

q i * x d x q 

(5) 

V <5 ‘ Vd 

(6) 

X„ = - i_x" 

q d q q 

(7) 

x„ (x -X") 

p _ q p»» + q q > 
e d x? e d + x" *q 

q q 

(8) 

_ x d’ x d . x d” x d „ 
e ql x^-xj e q " x^xj q 

(9) 

x d” x d . x d x d x d x d x d , 

e q2 ' xJj-xJ q xj xj-xj q xj d 

(10) 


The armature current Is given by 


* 
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The electromagnetic torque Is: 

T » * £ e SV e SV («) 

Equations (1-12) are In per unit based upon the machine rating. In these equations, 
the total tie line and transformer reactances. -X' , are Included In the machine 
Impedances. Therefore, voltages V d and V q In Eqs. (4) and (5) are the d-axls 
and q-axls components of the bus voltage. 

Equations of Motion 

The electromagnetic torque (referred to the rotor shaft) In lb-ft Is 

T L " T ^ (kv * k ««H738)T e (13) 

Two configurations for the equations of motion, or swing equations are 
formulated. The first representation of the drive train system Is by a ir-mechanical 
network (Fig. 5). 

The corresponding equations are: 


-ar ■ a - «n * a e 

(14) 

d4 



(15) 

<n_2K( V y-2T L -Bn 
dt J 

(16) 

dzw 

IF 3 

(17) 


where Is In mechanical radians and T Is In lb-ft. 

HI 

In the Hod-0, a Falk coupling Is Installed between the low speed shaft and 
the step up gearbox (not shown In Fig. 2). The torsional stiffness coefficient* 
Kf, of the Falk coupling Is a function of the power output* P Q » In kw. 

K f - (0.000042 P* + 0.0057 |P 0 | + 0.06) (18) 


w QVAun 


3 - 


The spring constant, K, In Fig. 5 is 

1 10 ® 

K r x ib-ft/rad . (19) 

1.919253 +tt- u 

*f 

If the Falk coupling is removed, 

K = 4.342 x 10 6 lb- ft/ rad . (20) 


The second configuration represents the drive train system by a ir-network 
with slip clutch, B r , Fig. 6. In this case, Eqs. (14) and (15) remain the same 
and Eqs. (16) and (17) are replaced by 


dn = T 2~ T L 
dt ” J 2 

m 2T-2K(« w -« r ) - B(n c +ZM) 
"dT " 


dZR 2K(W-2T C -B( V 2R) 
dt 


( 21 ) 

( 22 ) 

(23) 

(24) 


The Slip, S Is defined as 

S * 


fl N +0 e 

V ZR 


(25) 


where ^ and Q c are the reference shaft speeds (referred to the turbine rotor shaft) 
on the right hand and left hand sides of the slip clutch In Fig. 6, respectively. 

The mechanical torque transferred by the slip clutch In lb- ft is approximated 
by a function of the slip S: 


T 2 - as 4 + bS 3 + cS 2 + dS 

(26) 

T 2 

T c“T* 

(27) 


4 - 



Equations of the Blade Pitch Control System 

The Mod-0 feedsback the power output of the alternator to, the blade pitch 
control system to stabilize the output under load conditions. Here the pitch 
controller. Fig. 4, consists only of an Integrator, without the proportional 


components (K^ =0). 


The following equations are formulated according to Fig. 4. 


dX, 

TT ' MW 

(28) 

where 


P 0 - (kva base)(V d 1 d +V q 1 q ) 

(29) 

and P N Is the reference power input. 


n 

-< 

(30) 

dX, 

1 s X 

dt a 4 

(31) 

dX 4 2 2 

IT * “N X 2 " V3 ‘ 2 ^N X 4 

(32) 

The input and output variables of the flow limiter are: 


X 

II 

^1- 

X 

1 

(33) 

Y = limit ( - , 22 ,5 * x ) 

(31) 

The change In pitch angle Is 



(35) 

The pitch angle Is 


0 ■ + A0 

(35) 

The corresponding torque developed by the rotor In Ib-ft Is: 


T * f(v w ,n p ,0) 

(37) 
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Equation (37) is a nonlinear function of wind speed in mph, rotor speed 
in rpm, and pitch angle in mechanical degrees, [6]. The turbine blade pitch 
angle physically can vary only from -90 degrees to zero degrees so that the 
slope of the torque angle curves is always positive. Figure 7 shows the turbine 
torque-pitch angle curves for 15, 18 and 20 mph. 


Equations of Generator Excitation Control System 

A static excitation regulator is used to hold the terminal voltage, the 
output power factor, and the armature current of the generator within satisfactory 
limits. The excitation control system is approximated by the following equations: 


dt 


AV„ = 


dAV f 
“dt" = z 

(38) 

- = — - AV, - — 

T n t T 

e e 

(39) 

f = V fn + AV f 

(40) 

P. 

- K 3 (0.8 i a - v t ( kva baseP 

(41) 

td = V d " Vq 

(42) 

tq = V q + Vd 

(43) 

V t “7 V W +V Iq 

(44) 

■ s1n6, V q * cos6 

(45) 

tt 2 m 

(46) 


The first term in Eq. (41) corrects for the deviation of the generator 
terminal voltage, V t , from its reference value, V N - The second term corrects for 
power factor deviation. At rated output, the practical power factor is equal to 
0.8, lagging current. V t(J and In Eqs. (42) and (43) are the d-axis and q-axis 
components of V t . 
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DIGITAL SIMULATION 


The system equations formulated In the preceding sections are nonlinear 
and interrelated. They are solved numerically by the fourth order Adams-Bashforth 
(predictor) and Adams-Moulton (corrector) method, [7], for which a fixed stepslze 
of h *0.005 Is used. > 

Initially, the wind generator Is supplying rated power (100-kw) to the bus 
at 0.8 power factor, lagging current, and rated frequency. The tower shadow 
cuts In at time t=0. All case studies simulated here were run for 30 simulated 
seconds. Only the last 14 seconds of the simulation are plotted here. In order 
to allow the Initially stable system time to adjust to the artificially abrupt 
Imposition of the tower shadow at time*0. All time references on the graphs are 
relative to these last 14 seconds. 

System Input 

In all the cases under study, a constant wind speed of 18 mph Is assumed. 
However, practical experience with the Mod-0 has revealed that a significant 
cyclical perturbation occurs, caused by the obstruction of the wind by the 
supporting tower of r.n« turbine [8]. The rotor blades, which are downwind of 
the tower, enter Into this ‘tower shadow* of reduced wind velocity during each 
revolution. The shadow's arc Is estimated to be 30 degrees. 

Data from the NASA-Lewls Research Center reveals that the Mod-0 sustains 
approximately a 28 percent decrease In torque as a blade passes through the 
tower shadow. Fig. 8. The effective wind speed at the blades Is calculated by 
using Eq. (37) and holding both the rotor speed and the blade pitch angle constant. 
Thus the 28 percent decrease in torque Is equivalent to a 10 percent decrease In 
wind velocity at the blade when In the shadow. Fig. 9. A cosine function models 
the tower shadow effect. 
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Case Studies 



Case 1. Falk coupling* no slip clutch, Xl c 0.02 p.u. 

This case represents the Mod-0 as originally conceived, without a slip clutch. 
Early operational experience with this configuration was unsatisfactory. Later, 
a slip clutch was added to the Mod-0 (see Case 2). The configuration here 
Includes the Mod-O’s Falk coupling, and its existing tie-line and transformer 
reactances, X^. 

The following table summarizes the maximum and minimum values for the last 
14 simulated seconds of Case 1: 



Alternator 

Output 

(100 

Armature 

Current 

(p.u.) 

■Power 
Angle 
(El. Deg.) 

Rotor 

Speed 

(rpo) 

Pitch 

Angle 

(deg.) 

Turbine 

Torque 

(ft-lb) 

Generator 

Torque 

(ft-lb) 

Field 

yoltage 

(p.u.) 

Power 

Factor 

Maximum 

139.4 

1.305 

34.33 

40.12 

-0.267 

22898 

25295 

2.911 

0.852 

Minimum 

67.8 

0.733 

20.65 

B9 

EE2 

16420 

12167 

2.830 

0.715 

Max-Mi n 

71.6 

0.572 

13.68 

0.22 

0.257 

6478 

13128 

0.081 

0.137 
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Graphic representations for Case 1 are given In Figs. lOa-g. 

Unlike all the other cases, Case 1 was simulated for 40 seconds rather than 
30 seconds. Thus, as compared to the subsequent cases, an additional 10 seconds 
was provided here to allow this particularly unstable configuration to react to 
the abrupt Imposition of the tower shadow. 

Note that although the input Is periodically uniform, the graphs of this 
configuration reveal a moderate harmonic response. 

Case 2. Falk coupling, slip clutch, Xl= 0.02 p.u. 

This case represents the existing Mod-0 with its added slip clutch. 

The addition of the slip clutch has reduced the oscillations of all the 
variables tabulated below, except turbine torque, by a factor greater than 3. 
However, the response of the system here is characterized by a more severe 
harmonic response to the input. Figs. Ua-h. 

The following table summarizes the maximum and minimum values for the last 
14 simulated seconds of Case 2: 



Alternator 

Output 

<KW) 

Armature 

Current 

(p.u.) 

Power 
Angle 
(El. Deg.) 

Rotor 

Speed 

(rpm) 

Pitch 

Angle 

(deg.) 

Turbine 

Torque 

(ft-lb) 

Generator 

Torque 

(ft-lb) 

Field 

Voltage 

(p.u.) 

Power 

Factor 

Maxi mum 

112.2 

1.085 

29.70 

40.04 

+0.010 

23093 

20206 

2.878 

0.818 

Minimum 

89.5 

0.906 

25.47 

39.97 

-0.063 

16632 

16076 

2.856 

0.779 

Max-Mi n 

22.7 

0.179 

4.23 

0.07 

0.073 

6461 

4130 

0.022 

0.039 


Case 3. No Falk coupling, no slip clutch, XfQ.02 p.u. 

This case is a theoretical one, differing from Case 1 only in the removal of 
the Mod-O's Falk coupling. The improvement in performance caused by this 
hypothetical stiffer shaft is considerable. The difference between maximum and 
minimum values here is very similar in magnitude to the corresponding values of 
Case 2. However, contrary to the harmonic response exhibited in the simulation 
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of the Mod-0 with a slip clutch (Case 2), maximum and minimum values for each 
cycle in Case 3 are very nearly constant. Figs. 12a-g. 

The following table summarizes the maximum and minimum values for the 
last 14 simulated seconds of Case 3: 



A1 temator 
Output 
0W) 

Armature 

Current 

(p.u.) 

Power 
Angle 
(El. Deg.) 

Rotor 

Speed 

(rpm) 

Pitch 

Angle 

(deg.) 

Turbine 

Torque 

(ft-lb) 

Generator 

Torque 

(ft-lb) 

field 

Voltage 

(p.u.) 

Power 

Factor 

Maximum 

108.4 

1.055 

29.34 

40.04 

-0.367 

22759 

19533 

2.880 

0.820 

Minimum 

85.4 

0.871 

24.78 

39.95 

-0.448 

16449 

15352 

2.857 

0.775 

Max-Mi n 
— 

23.0 

0.184 

• 4.56 

0.09 

0.081 

6310 

4181 

0.023 

.045 


Case 4. No Falk coupling, slip clutch, X|_ =0,02 p.u. 

This hypothetical case attempts to combine the successful methods of 
Cases 2 and 3. The configuration here differs from the existing Mod-0 only 
in the removal of the Falk coupling. Differences between maxima and minima for 
the variables tabulated below (except for wind turbine torque) are approximately 
one half of the corresponding values for Case 2 or Case 3. Additionally, the 
harmonics which appeared in Case 2 have been eliminated by removing the Falk 
coupling. 

Figure 13a-h are the plots obtained from Case 4. 

The following table summarizes the maximum and minimum values for the 
last 14 simulated seconds of Case 4: 



A1 temator 
Output 
(KM) 

Armature 

Current 

(p.u.) 

Power 
Angle 
(El. Deg.) 

Rotor 

Speed 

(rpm) 

Pitch 

Angle 

(deg.) 

Turbine 

Torque 

(ft-lb) 

Generator 

Torque 

(ft-lb) 

B 

Power 

Factor 

Maximum 

105.7 

1.036 

28.48 

40.02 

-0.006 

23075 

19021 

2.873 

0.807 

Minimum 

93.4 

0.937 

26.26 

39.98 

-0.046 

16640 

16790 

2.861 

0.788 

Max-Mi n 

— - -- 

12.3 

.099 

2.22 

0.04 

0.040 

6435 

2231 

0.012 

0.019 
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The subsequent Cases 5-9 attempt to refine the hypothetical Case 4, i.e.. 
Cases 5-9 simulate the Mod-0 with the slip clutch and without the Falk coupling. 
However, In each of these cases the tie line and transformer reactances, X^, has 
a different value. 

The previous simulations revealed that the terminal voltage, V t , was nearly 
equal to the constant bus voltage. Under such conditions, the field excitation 
voltage, V f , cannot be effectively controlled by the voltage regulator when 
necessary. Since the power factor, armature current, and power angle must be 
held within narrow limits, the patent means of adjusting to reduce output 
oscillations caused by the tower shadow is to vary within an acceptable range. 
Changing X^ will not only allow V^. to vary within reasonable limits, but also 
can store or release energy to smooth the output oscillations. 

Cases 4-9 simulate the Mod-0 with only the value of X^ (and hence the initial 
conditions as well) different in each case. 

Case 5. No Falk coupling, slip clutch, Xi = 0.40 p.u. 

As the plots (Fig. 14a, b) show, the deviations from rated alternator output 
have been reduced by about 20 percent, as compared to Case 4, as a result of 
Increasing X^ from 0.02 to 0.40 p.u. The differences between maximum and minimum 
values of armature current, blade pitch angle, generator torque and power factor 
have also been reduced. However, the field voltage and rotor speed have both 
Increased their differences between maximum and minimum values. The former is 
not surprising since the purpose of increasing X^ is improved performance of the 
static excitation control system. 

The power angle has Increased the magnitude of the difference between its 
maximum and minimum values by nearly 50 percent. However, the Increase in X L 
has also resulted In a significant increase in the mean value of the power angle. 
Therefore, let us consider the relative change in the power angle, 
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n/«/ _ maximum (power angle) - minimum (power angle) 

RC( power angle) - ( pwef arTgl elT/2 


We find that RC(power angle) of Case 4 = 0.0811 electrical degrees, and 
RC(power angle) of Case 5 =0.0787 electrical degrees. Using the relative 
change criterion, the power angle has less fluctuation here than in Case 4. 
The relative change criterion for all other variables is consistent with the 
difference between maximum and minimum values. 


The following table summarizes the maximum and minimum values for the last 
14 simulated seconds of Case 5: 



Alternator 

Output 

(KW) 

Armature 

Current 

(p.u.) 

Power 
Angle 
(El. Deg.) 

Rotor 

Speed 

(rpm) 

Pitch 

Angle 

(deg.) 

Turbine 

Torque 

(ft-lb) 

Generator 

Torque 

(ft-lb) 

Field 

Voltage 

(p.u.) 

Power 
Factor 


Maximum 

104.4 

0.888 

39.88 

40.03 

-0.067 

22968 

18693 

2.723 

E3 

Minimum 

94.7 

0.821 

36.86 

39.98 

- 0.100 

16571 

16938 

2.710 

0.791 

Max-Mi n 

9.7 

0.067 

3.02 

0.05 

0.033 

6397 

1755 

0.013 

0.015 


Case 6. No Falk coupling, slip clutch, Xt=0.47 p.u. 

In this case, is again increased. As the plots of Case 6 (Fig. 15a-h) 
show, maximum and minimum values for each cycle are still reasonably constant, 
although a very slight harmonic can be detected. As compared to Case 5, perfor- 
mance has improved here with regard to the alternator output, the armature 
current, the blade pitch angle, the generator torque, the power factor, and the 
power angle (using the relative change criterion, RC(power angle) of Case 6 is 
0.0776 electrical degrees). Performance has been degraded marginally with respect 
to the rotor speed, and the field voltage. 

The following table summarizes the maximum and minimum values for the last 
14 simulated seconds of Case 6: 


ORIGINAL 
OF POOR QUALITY 
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Alternator 

Output 

(XW) 

Armature 

Current 

(p.u.) 

Power 
Angle 
(El. Deg.) 

Rotor 

Speed 

(rpm) 

Pitch 

Angle 

(deg.) 

Turbine 

Torque 

(ft-lb) 

Generator 

Torque 

(ft-lb) 

Field 

Voltage 

(p.u.) 

Power 

Factor 

Maximum 

104.2 

0.872 

42.06 

40.03 

-0.073 

22958 

18640 

2.715 

0.805 

Minimum 

95.0 

HB 

38.92 

39.98 

-0.105 

16563 

16981 

2.693 

0.791 

Max-Mi n 

9.2 

0.063 

3.14 

0.05 

0.032 

6395 

1659 

0.022 

0.014 


Case 7. No Falk coupling, slip clu tch, X =0.50 p.u. 

In this case, the small increase in X L has resulted in marginal improvements 
over Case 6 with respect to the alternator output, armature current, blade pitch 
angle, and generator torque. However, performance fell off somewhat with respect 
to the power angle (using either absolute difference or relative change criteria), 
the rotor speed, the field voltage and the power factor. Furthermore, as the 
plots of Case 7 (Fig. 16a, b) show, the slight harmonic observed in Case 6 persists 
in Case 7. 


The following table summarizes the maximum and minimum values for the last 
14 simulated seconds of Case 7: 



A1 temator 
Output 
(KW) 

Armature 

Current 

(p.u.) 

Power 
Angle 
(El . Deg.) 

Rotor 

Speed 

(rpm) 

Pi tch 
Angle 
(deg.) 

.Turbine 

Torque 

(ft-lb) 

Generator 

Torque 

(ft-lb) 

Field 

Voltage 

(p.u.) 

Power 

Factor 

Maximum 

104.1 

0.867 

43.03 

40.03 

-0.076 

22954 

18619 

2.712 

0.806 

Minimum 

95.0 

0.805 

39.79 

39.97 

-0.107 

16561 

16984 

2.688 

0.791 

Max-Mi n 

9.1 

0.062 

3.24 

0.06 

0.031 

6393 

1635 

0.024 

0.015 


Case 8. No Falk coupling, slip clutch, X|_=0.55 p.u. 

In this case, the Increase in X L has resulted in smaller differences between 
maximum and minimum values of alternator output and generator torque as compared 
to the corresponding values of Case 7. However, the differences between maximum 
and minimum values of field voltage, rotor speed and power angle have continued 
to Increase over the previous case. Also note In the plots (Fig. 17a, b) that the 
harmonics have become more conspicuous. 
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The following table summarizes the maximum and minimum values for the last 
14 simulated seconds of Case 8: 



Alternator 

Output 

m 

Armature 

Current 

(P.u.) 

Power 
Angle 
(El. Deg.) 

Rotor 

Speed 

(rpm) 

Pitch 

Angle 

(deg.) 

Turbine 

Torque 

(ft-lb) 

Generator 

Torque 

(ft-lb) 

Field 

Voltage 

(p.u.) 

Power 

Factor 

Maximum 

104.0 

0.858 

44.71 

40.03 

-0.079 

22949 

18591 

2.715 

0.806 

Minimum 

95.1 

0.797 

41.27 

39.97 

-0.110 

16557 

16994 

2.672 

0.791 

Max- Min 

8.9 

0.061 

3.44 

0.06 

0.031 

6392 

1597 

0.043 

0.015 


Case 9. No Falk coupling, sli p clutch, Xi s 0.60 p.u. 

As the plots of Case 9 (.-ig. 18a, b) s^ow, the harmonics have here become 
quite prominent. Therefore, simulation of further Increases of with this 
configuration was abandoned after Case 9. 

As compared to Case 8, the increase in here did result with a smaller 
difference between maximum and minimum values of alternator output, armature 
current, and generator torque. However, the differences were larger with respect 
to the power angle, rotor speed, pitch angle, field voltage and power factor. 

The following table summarizes the maximum and minimum values for the last 
14 simulated seconds of Case 9: 



Alternator 

Output 

W 

Armature 

Current 

(p.u.) 

Power 
Angle 
(El. Deg.) 

Rotor 

Speed 

(rpm) 

Pitch 

Angle 

(deg.) 

Turb.ne 

Torque 

(ft-lb) 

Generator 

Torque 

(ft-lb) 

Field 

Voltage 

(p.u.) 

Power 

Factor 

Maximum 

104.0 

0.850 

46.47 

40.03 

-0.081 

22946 

18574 

2.726 

0.807 

Minimum 

95.2 

0.791 

42.75 

39.97 

-0.113 

16553 

16998 

2.650 

0.791 

Max-Mi n 

8.8 

0.059 

3.72 

0.06 

0.032 

6393 

1576 

0.076 

0.016 


As Case 5-9 demonstrate, changes to the tie-line and limiting reactances, 
X^, can significantly Improve the overall performance of the hypothetical 
Mpd-0 (l.e., the Mod-0 without a Falk coupling). How would such a change In 
X^ effect the existing Mod-0? 


Case 10. Falk coupling, slip clutch, X|_ * 0.47 p.u. 


A comparison of the plots obtained In this case. Fig. 19a-h, with those of 
the simulation of the existing Mod-0 (Case 2), reveals no significant enhancement 
of the stabilization process. For the 2 cases, differences between maximum and 
minimum values of alternator output, blade pitch angle, and generator torque, 
as well as armature current and power angle (relative change), are approximately 
the same. The rotor speed for Case 10 Is significantly worse. Moreover, the 
harmonics which plagued Case 2 have become further aggravated with the Increase 
in X^. 

The following table summarizes the maximum and minimum values for the last 
14 simulated seconds of Case 10: 



Alternator 

Output 

(KM) 

Armature 

Current 

(p.u.) 

Power 
Angle 
(El. Deg.) 

Rotor 

Speed 

(rpm) 

Pi tch 
Angle 
(deg.) 

Turbine 

Torque 

(ft-lb) 

Generator 

Torque 

(ft-lb) 

Field 

Voltage 

(p.u.) 

Power 

Factor 

Maximum 

111.4 

0.922 

44.46 

40.07 

-0.054 

22965 

19940 

2.723 

0.814 

Minimum 

89.1 

0.769 

36.93 

39.94 

-0.125 

16561 

15917 

2.690 

0.782 

Max-Mi n 

22.3 

0.153 

7.53 

0.13 

0.071 

6405 

4023 

0.033 

0.032 


The remaining Case 11-18 use the previous technique of parameterization 
of X^, using Case 3 as a basis. Recall that Case 3 achieved substantial success 
by simply removing the Falk coupling and slip clutch. 

Case 11. No Falk coupling, no slip clutch, X( > s 0.20 p.u. 

The plots of Case 11, Fig. 20a, b, clearly demonstrate the Increased 
stability brought about by the 10 fold Increase In X^. In Case 11, differences 
between maximum and minimum values for all variables listed (except turbine 
torque) In the table below are about one third less than the corresponding values 
In Case 3. Moreover, maximum and minimum values for Case 11 remain very nearly 
constant In each cycle. 
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The following table summarizes the maximum and minimum values for the last 


14 simulated seconds of Case 11: 



Alternator 

Output 

(KW) 

Armature 

Current 

(p.u.) 

Power 
Angle 
(El. Deg.) 

Rotor 

Speed 

(rpm) 

Pitch 
Angle 
(deg. ) 

Turbine 

Torque 

(ft-lb) 

Generator 

Torque 

(ft-lb) 

Field 

Voltage 

(p.u.) 

Power 

Factor 

Maximum 

105.6 

0.953 

34.18 

40.04 

-0.408 

22681 

18963 

2.779 

0.812 

Minimum 

89.9 

0.839 

30.17 

39.96 

-0.467 

16408 

16111 

2.766 

0.781 

Max-Mi n 

15.7 

0.114 

br 

0.08 

0.059 

6273 

2852 

0.013 

0.031 


Case 12. No Falk coupling, no slip clutch, Xl = 0.30 p.u. 

The plots. Fig. 21a, b, and the table below illustrate the enhanced stability 
of the alternator output, armature current, power angle, and generator torque. 

The peak values of the plots remain very nearly constant. 

The following table summarizes the maximum and minimum values for the 
last 14 simulated seconds of Case 12: 



Alternator 

Output 

(KM) 

Armature 

Current 

(p.u.) 

Power 
Angle 
(El. Deg.) 

Rotor 

Speed 

(rpm) 

Pitch 

Angle 

(deg.) 

Turbine 

Torque 

(ft-lb) 

Generator 

Torque 

(ft-lb) 

Field 

Voltage 

(p.u.) 

Power 

Factor 

Maximum 

104.9 

0.933 

36.03 

40.04 

-0.416 

22666 

18833 

2.822 

0.787 

Minimum 

91.2 

0.840 

32.16 

39.96 

-0.468 

16402 

16351 

2.808 

0.758 

Max- Min 

13.7 

0.093 

3.87 

0.08 

0.052 

6264 

2482 

0.014 

0.029 


Case 13. No Falk coupling, no slip clutch, Xi*0.35 p.u. 

Plots of Case 13, Fig. 22a, b, demonstrate the further improvement of 
system response. 

The following table summarizes the maximum and minimum values for the 
last 14 simulated seconds of Case 13: 
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Alternator 

Output 

(KM) 

Armature 

Current 

(p.u.) 

Power 
Angle 
(El. Deg.) 

eJjgH 

Pitch 

Angle 

(deg.) 

Turbine 

Torque 

(ft-lb) 

Generator 

Torque 

(ft-lb) 

Bi 

Power 

Factor 

Maximum 

104.7 

0.924 

37.05 

40.04 

-0.420 

22660 

18789 

2.844 

0.777 

Minimum 

91.7 

0.838 

33.14 

39.96 

-0.469 

16397 

16438 

2.821 

0.748 

Max-Mi n 

13.0 

0.086 

3.91 

0.08 

0.049 

6263 

2351 

0.023 

0.029 


Case 14. No Falk coupling, no slip clutch, Xl= 0.40 p.u. 

Plots of Case 14, Fig. 23a, b, illustrate the improvement of system response. 
The following table summarizes the maximum and minimum values for the 
last 14 simulated seconds of Case 14: 



Alternator 

Output 

(KW) 

Armature 

Current 

(p.u.) 

Power 
Angle 
(El. Deg.) 

Rotor 

Speed 

(rpm) 

Pitch 

Angle 

(deg.) 

Turbine 

Torque 

(ft-lb) 

Generator 

Torque 

(ft-lb) 

Field 

Voltage 

(p.u.) 

Power 

Factor 

Maximum 

104.4 

0.888 

39.95 

40.04 

-0.434 

22631 

18705 

2.730 

0.808 

Minimum 

92.4 

0.806 

36.14 

39.96 

-0.480 

16379 

16541 

2.705 

0.787 

Max-Mi n 

12.0 

0.082 

3.81 

0.08 

0.046 

6252 

2164 

0.025 

0.021 


Case 15. No Falk coupling, no slip clutch, Xl= 0.43 p.u. 

Plots of Case 15, Fig. 24a-g, demonstrate the continued reduction of the 
difference between extreme values of all variables, except rotor speed and 
field voltage. System response here is comparable to the case with a slip clutch, 
no Falk coupling, and with still 0.02 p.u. (Case 4). However, a very slight 
harmonic response is noticeable here, particularly at the low end of the time 
scale. 

The following table summarizes the maximum and minimum values for the 
last 14 simulated seconds of Case 15: 


ORIGINAL 
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Alternator 

Output 

(KW) 

Armature 

Current 

(p.u.) 

Power 
Angle 
(El. Oeg.) 

Rotor 

Speed 

(rpm) 

Pitch 

Angle 

(deg.) 

Turbine 

Torque 

(ft-lb) 

Generator 

Torque 

(ft-lb) 

Field 

Voltage 

(p.u.) 

Power 

Factor 

Maximum 

104.3 

0.881 

40.89 

40.04 

-0.437 

22627 

18681 

2.731 

0.808 

Minimum 

92.7 

0.801 

37.04 

39.96 

-0.483 

16375 

16586 

2.695 

1223 

Max-Mi n 

11.6 

0.080 

3.85 

0.08 

0.046 

6252 

2095 

0.036 

0.021 


Case 16. No Falk coupling, no slip clutch, Xl= 0.47 p.u. 

In this case, X L Is Increased by only 0.04 p.u. over the previous Case 15. 

The plots of Case 16, Fig. 25a-g, reveal a slight reduction In the differences 
between maximum and minimum values of alternator output, armature current and 
generator torque from the corresponding values in Case 15. However, the power 
angle, power factor, and field voltage were better in Case 15. Host significantly, 
the slight harmonics which appeared at the low end of the time scale In Case 15 
extend over the entire time scale In Case 16, and are of greater amplitude. 

The following table surnnarizes the maximum and minimum values for the last 
14 simulated seconds of Case 16: 



Alternator 

Output 

(KW) 

Armature 

Current 

(p.u.) 

Power 
Angle 
(El. Deg.) 

Rotor 

Speed 

(rpm) 

Pitch 

Angle 

(deg.) 

Turbine 

Torque 

(ft-lb) 

Generator 

Torque 

(ft-lb) 

Field 

Voltage 

(p.u.) 

Power 

Factor 

Maximum 

104.3 

0.873 

42.27 

40.04 

-0.439 

22622 

18678 

2.739 

0.809 

Minimum 

92.8 

0.795 

38.12 

39.96 

-0.486 

16369 

16609 

2.677 

0.786 

Max-Mi n 

11.5 

0.078 

4.15 

0.08 

0.047 

6253 

2069 

0.062 

0.023 


Case 17. No Falk coupling, no slip clutch, Xi»0.50 p.u. 

Plots of Case 17, Fig. 26a,b, show that the harmonic response of the 
system has become aggravated with this most recent increase In X^. 

The following table surnnarizes the maximum and minimum values for the 
last 14 simulated seconds of Case 17: ' 
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Alternator 

Output 

(KW) 

Armature 
Current 
(p.u ) 

Power 
Angle 
(El. Deg.) 

Rotor 

Speed 

(rpm) 

Pitch 

Angle 

(deg.) 

Turbine 

Torque 

(ft-lb) 

Generator 

Torque 

(ft-lb) 

Field 

Voltage 

(p.u.) 

Power 

Factor 

Maximum 

104.4 

0.867 

43.46 

40.04 

-0.439 

22623 

18697 

2.754 

0.812 

Minimum 

92.8 

0.791 

38.81 

39.96 

-0.490 

16364 

16610 

2.649 

0.784 

Max-Mi n 

11.6 

0.076 

4.65 

0.08 

0.051 

6259 

2087 

0.105 

0.028 


Case 18. No Falk coupling, no slip clutch, Xi=0.55 p.u. 

Plots of Case 18, Fig. 27a, b, clearly demonstrate the degenerate response 
of the control system a* caused by the Increased value of X L * Although the 
system is still functioning within acceptable bounds, the difference between 
extreme values for all tabulated variables has increased to about the same 
level as Case 12 (5^=0. 20). Unlike Case 12, however, the extreme values for 
each cycle oscillate significantly, no longer reflecting the periodic input. 

The difference between maximum and minimum values of the field voltage is 
greater in Case 18 than in any of the previous cases, while the power factor 
is only surpassed in Case 1. 

The following table summarizes the maximum and minimum values for the last 
14 simulated seconds of Case 18: 



Alternator 

Output 

(KM) 

Armature 

Current 

(p.u.) 

Power 
Angle 
(El. Deg.) 

Rotor 

Speed 

(rpm) 

Pitch 

Angle 

(deg.) 

Turbine 

Torque 

(ft-lb) 

Generator 

Torque 

(ft-lb) 

Field 

Voltage 

(p.u.) 

Power 

Factor 

Maximum 

105.5 

0.862 

46.21 

40.05 

-0.427 

22642 

18876 

2.868 

0.823 

Minimum 

91.9 

0.783 

39.33 

39.95 

-0.506 

16342 

16435 

2.522 

0.771 

Max-Mi n 

13.6 

0.079 

6.88 

0.10 

0.079 

6300 

2441 

0.346 

0.052 
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Case 

Falk 

Coupling 

Slip 
Cl utch 

*L 

(p.u.) 

Alternator Output 
(KW) 

Armature Current 
(p.u.) 

Power Angle 
(el. deg.) 

Rotor Speed 
(rpm) 

Generator Torque 
(ft-lb) 

max 

min 

max- 
mi n 


— 

min 

max- 
mi n 

max 

min 

max- 
mi n 

max 

min 

max- 
mi n 

max 

min 

max- 

min 

1 

Yes 

No 

0.02 

B 

67.8 

71.6 

1.305 

0.733 

0.572 

34.33 

20.65 

13.68 

40.12 

39.90 

0.22 

25295 

12167 

13128 


Yes 

Yes 



89.5 

22.7 

1.085 

0.906 

0.179 

29.70 

25.47 

4.23 

40.04 

39.97 

0.07 

20206 

16076 

4130 


No 

No 



85.4 

23.0 

1.055 

0.871 

0.184 

29.34 

24.78 

4.56 

40.04 

39.95 

0.09 

19533 

15352 

4181 


No 

*Yes 


105.7 

93.4 

12.3 

1.036 

0.937 

0.099 

28.48 

26.26 

2.22 

40.02 

39.98 

0.04 

19021 

16790 

2231 


No 

Yes 


104.4 

94.7 

9.7 

0.888 

0.821 

0.067 

39.88 

36.86 

3.02 

40.03 

39.98 

0.05 

18693 

16938 

1755 


No 

Yes 


104.2 

95.0 

9.2 

0.872 

0.809 

0.063 

42.06 

3b. 92 

3.14 

40.03 

39.98 

0.05 

18640 

16981 

1659 


No 

Yes 

BBS 

104.1 

95.0 

9.1 

0.867 

0.805 

0.062 

43.03 

39.79 

3.24 

40.03 

39.97 

0.06 

18619 

16984 

1635 


No 

Yes 

; 

104.0 

95.1 

8.9 

0.858 

0.797 

0.061 

44.71 

41 .27 

3.44 

40.03 

39.97 

0.06 

18591 

16994 

1597 


No 

Yes 

0.60 

104.0 

95.2 

8.8 

0.850 

0.791 

0.059 

46.47 

42.75 

3.72 

40.03 

39.97 

0.06 

18574 

16998 

1576 

10 

Yes 

Yes 

0.47 

111.4 

89.1 

22.3 

0.922 

0.769 

0.153 

44.46 

36.93 

7.53 

40.07 

39.94 

0.13 

19940 

15917 

4027 

11 

No 

No 

0.20 

105.6 

89.9 

15.7 

0.953 

0.839 

0.114 

34.18 

30.17 

4.01 

40.04 

39.96 

0.08 

18963 

16111 

2852 

12 

No 

No 

0.30 

104.9 

91.2 

13.7 

0.933 

0.840 

0.093 

36.03 

32.16 

3.87 

40.04 

39.96 

0.08 

18833 

16351 

2482 

13 

No 

No 

0.35 

104.7 

91.7 

13.0 

0.924 

0.838 

0.086 

37.05 

33.14 

3.91 

40.04 

39.96 

0.08 

18789 

16438 

2351 

14 

No 

No 

0.40 

104.4 

92.4 

12.0 

0.883 

0.806 

0.082 

39.95 

36.14 

3.81 

40.04 

39.96 

0.08 

18705 

16541 

2164 

IS 

No 

No 

0.43 

104.3 

92.7 

11.6 

0.881 

0.801 

0.080 

40.89 

37.04 

3.85 

40.04 

39.96 

0.08 

18681 

16586 

2095 

16 

No 

No 

0.47 

104.3 

92.8 

11.5 

0.873 

0.795 

0.078 

42.27 

33.12 

4.15 

40.04 

39.96 

0.08 

186 78 


2069 

17 

No 

No 

0.50 

104.4 

92.8 

11.6 

0.867 

0.791 

0.076 

43.46 

38.81 

4.65 

40.04 

39.96 

0.08 

18697 

il 

2087 

18 

No 

No 

0.55 

105.5 

91.9 

13.6 

0.862 

0.783 

0.079 

46.21 

39.33 

6.88 

40.05 

39.95 

0.10 

18876 

16435 

2441 


Note: Case 1 was simulated for 40 seconds. All other cases were simulated for 30 seconds. 


MAXIMUM AND MINIMUM VALUES FOR CASES 1-18 































CONCLUSIONS 


1. Removal of the Falk coupling results in a stiffer shaft, which 
dramatically increases system stability. Its removal is therefore 
recommended. 

2. The slip clutch effectively aids system stability both with and withou. 
the Falk coupling. Once the Falk coupling is removed, output oscillations 
caused by the tower shadow (28 percent decrease in turbine torque) can 

be held by the slip clutch to within 7 percent of the rated value. 

However, the slip clutch is an expensive method of control. Furthermore, 
the clutch makes it impractical for implementation by the larger machines 
yet to be built by ERDA-NASA. 

3. The reactances X L significantly affects the dynamic stability. For values 
of the tie line and transformer reactance, X L , at the level currently 
being used in the Mod-0, the excitation control system has little effect 
on stabilization. Additional reactance (or a longer transmission line) 
is needed to permit the static excitation control system to realize its 
control capabilities. The simulations reveal that increasing X^ to 
within a narrow range around 0.43 p.u. (for the configuration without a 

a slip clutch), or 0.47 p.u. (for the configuration with a slip clutch), 
can reduce the variation of power output caused by the tower shadow 
(for a hypothetical Mod-0 without a Falk coupling) from 12 percent (Case 2) 
for the current Mod-0 to 7 percent (without a slip clutch) or 5 percent 
(with a slip clutch) of the rated value. 

Optimization of the tie line and transformer reactances is an 
inexpensive, easily implemented technique which can be applied to both 
large and small machines. 


4. As Case 10 demonstrates, the technique of using optimal does not 
appear to be compatible with a drive train configuration containing a 
Falk coupling. For that reason, it is recommended that not be 
increased to the range discussed above for the existing Mod-0, which has 
a Falk coupling. 

5. For values of X L less than the optimal range, control system effectiveness 
is reduced. However, the system does reflect the constant amplitude of 
the periodic input. 

For values of X^ exceeding the optimal range, the system no longer 
reflects the constant amplitude of the periodic input, becoming beset 
with severe harmonics. 

6. Comparison of the effects of the slip clutch (Case 4) with that of 
optimization of X^ (Case 15) reveals very similar variations of power 
output for the two techniques. Therefore, by virtue of its lesser 
expense, greater durability and wider range of applications, the optimi- 
zation of the tie-line and transformer reactance seems the better method 
for smoothing the output power of a wind turbine generator. 

7. Combining the optimal value of X^ with a slip clutch yields the best 
method for stabilizing the system (Case 6). However the small improve- 
ments realized by the installation of the slip clutch (Case 6 vs. 

Case 15) do not seem to warrant the added expense of the fluid drive. 

It is therefore recommended that the technique advanced here -- 
a static excitation control system with an optimal value of X L 
and removal of the Falk coupling — replace the slip clutch in 
future wind turbine generators. 
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NOMENCLATURE 


The symbols used correspond insofar as possible to standard definitions. 
Other symbols used are 

Xf.Xnq flux linkage of field and damper circuit in j-axis 

M f ,M maximum mutual Inductance between armature circuit and field {or 
J damper circuit in j axis) 

Lff»l_ni self Inductance of field circuit and damper circuit in j axis, 

J respectively 

a) base angular velocity 

P number of poles of the synchronous machine 

N gear ratio between high and low speed shaft 

Vf n ,V f reference or instantaneous field voltage referred to the stator 

B damping coefficient of the complete rotating system in lb-ft-sec 

J inertia of the complete rotating system (referred to low speed shaft) 

in lb-ft-sec2 

ft angular speed of alternator shaft (referred to wind turbine rotor 

shaft) in radians per second 

6 . angular displacement of j section of the rotor shaft in mechanical 

J radians from selected reference axis 

controller proportional gain in seconds 

l <2 controller integral gain per kw-sec 

9 pitch angle of blade in mechanical radians 


0 N reference pitch angle of blade in mechanical radians 

£ damping ratio 


ip control hydraulic actuator time constant in seconds 

a> N undamped natural angular frequency in radians per second 
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SYSTEM CONSTANTS 


R a =0.018, x^'2.21, x^=1.064, x^ =0.165 p.u. , 

x^ = 0.128, x|| = 0.193, V b = 1.000 p.u., 

T^ 0 = 1 .94212, T^ 0 = 0.01096, T" 0 = 0.06230 sec 

N = 45, B =959 Ib-ft-sec 

J = 101,100, =95,760, J 2 =5340 lb-ft-sec 2 

Ki=0, < 2 = 6 x 10" 4 per kw-sec 

w N = 100 radians per sec, sec, £=0.02 

M e = 30.000, T e =0.0500 sec, <2 = 0.2000 

a = 7.1x!0 10 , b=3.3xl0 9 , c = -l.lxl0 8 , d = 1.8x!0 6 

f2.. =4.189, 51 =4.257 radians per sec 
N c 

AC<NOWLEDGMENTS 

We are grateful to Dr. W.A. Lewis for his valuable suggestions. This 
work is supported by NASA-Lewis Research Center and the Hawaii Natural 
Energy Institute. Thanks to the Center for Engineering Research for typing 
this manuscript. 


- 24 - 


ORIGINAL PAGE IS 
OF POOR QUALITY 



REFERENCES 


1. Ronald L. Thomas, "Large Experimental Wind Turbines— Where We are Now," 

NASA TMX-71 890. 

2. Craig C. Johnson and Richard T. Smith, "Dynamics of Wind Generators on 
Electric Utility Networks," IEEE Trans, on Aerospace and Electronic Systems, 

Vol . AES 12, No. 4, July 1976, pp. 483-493. 

3. R.H. Park, "Two Reaction Theory of Synchronous Machines, Generalized Method 
of Analysis, Part 1," AIEE Transactions , Vol. 48, July 1929, pp. 716-730. 

4. David W. Olive, "Digital Simulation of Synchronous Machine Transients," 

IEEE Trans, on Power Apparatus and Systems, Vol. PAS-87, Auqust 1968, 
ppT 1669-1975. 

5. H.H. Hwang and Leonard J. Gilbert, "Synchronization of Wind Turbine Generators 
Against an Infinite Bus Under Gusting Wind Conditions," presented at IEEE Power 
Engineering Society 1977 Summer Meeting and to appear in IEEE Trans, on Power 
Apparatus and Systems , paper no. F77675-2. 

6. Robert E. Wilson and Peter B.S. Lissman, Applied Aerodynamics of Wind Power 
Machines , Oregon State University, May 1975. 

7. Germund Dahlquist, Ake Brock, and translated by Ned Anderson, Numeri cal 

Methods , New Jersey: Prentice-Hall, 1974, Chapter 8, pp. 347-348. 

8. John C. Glasgow and Bradford S. Linscott, "Early Operation Experience on 
ERDA/NASA 100-kw Wind Turbine," NASA TMX-71 601. 


- 25 - 


APPENDIX: DERIVATIONS OF INITIAL CONDITIONS 


Synchronous Machines 

The following conditions can be found in Fig. 28. 


<}>b = } C 4> t + sin"’(sin<J> t cos<J> t )] 

v b 


■1 


2P_X 


o L 


I* = 


a V b cosd» b 
. v b cos *b 


v t = 


cosd> . 


. . -l(f ( V X L )cos *b ' R a sin +b ]I a l 

6 = U " K 7 [( V X L ,Sln '*b * 

V f = V b cos6 + R a I a (6+<b b ) + (X d +X L )I a sin'6+4. b ) 

’d = ! a Sln ( s+ *b> 

! q = I a cos(«+* b ) 


(47) 

(48) 

(49) 

(50) 

(51) 

(52) 
(51) 


Under steady-state, the derivatives of Eqs. (1), (2), and (3) are zero, and 


initial values of e d> e q2 » and e^, are 


e d * e q2 * 0 

*ql * V f 


Solving Eqs. (6-10) simultaneously, there yields 

*;--<v x d»d + v f 

•d* 


X .-XL X\-X*i 

p' s (\l + — — p») _iL 

e q l VxJ-xJ V x^-x” 


(52) 

(53) 

(54) 

(55) 

(56) 


d d 


Mechanical Networks 


6 W is found in Fig. 5 


6- = 


it 6 


m " 90" NF 


60, 


5 = V + - 

w K m 


N 

2K 


(57) 

(58) 
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( 59 ) 


The following initial values are easily derived from Fig. 6 

T - T c + ® c 

s r = s m 

Bn 

{ w-7 +s r'ir 


(60) 

(61) 
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Fig. 2 100-kw Wind Turbine Drive Train 
Assembly and Yaw System 



- 29 - 





I 

CO 


o 


I 


PITCH 



CONTROLLED 

VARIABLE 


WIND 

SPEED 


Fig. 4 The Pitch Control System of ERDA-NASA 100-kw Wind Turbine 
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Fig. 6 IT-mechanical Network with Slip Clutch 
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Fig. 7 The Turbine Torque-Pitch Angle Curve 
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Fig. 9 Wind Speed vs. Time, 

corresponding to Fig. 8 
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DUE. FT-LBS -• . WIND TURBINE TORQUE. FT-LBS 














Fig. 13d: Case 4, Rotor Speed 

No Falk coupling, slip clutch, X L = 0.02 p.u. 
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ROTOR SLIP 

0.01560 0.01592 0.01624 0.01656 0.01688 3 16700 




ROTOR SPEED, RPM 




Fig. 15a: Case 6, Alternator Output 

No Falk coupling, slip clutch, X L =0.47 p.u. 
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ROTOR SPEED, RPM 



r 



ROTOR SLIP 











Fig. 19c: Case 10, Generator Pov 

• Falk coupling, slip cl 



Fig. 19d: Case 10, Rotor Speed 

Falk coupling, slip c 
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GENERATOR POWER ANGLE, ELEC. DEGREES 







Fig. 24g: Case 15, Generator Torque 

No Falk coupling, no slip clutch, X L =0.43 p.u. 
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ARMATURE CURRENT. P.U 





WIND TURBINE BLADE PITCH ANGLE, DEGREES 
TORQUE, FT-LBS 2 _ n 490 -0.474 -0.458 -0.442 -0.426 



. 25e: Case 16, Blade Pitch Angle 

No Falk coupling, no slip clutch, Xj_ = 0.47 p.u. 
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Fig. 25f : Case 16, Turbine Torque 

No Falk coupling, no slip clutch, X L = 0.47 p.u. 
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Fig. 28 Phasor Diagram of Synchronous Generator 
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